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O Abstract 

Evaluating two-loop integrals in the transverse momentum space we ob- 
tain the trajectory of the Reggeized gluon in QCD in an explicit form in the 
two-loop approximation. It is presented as an expansion in powers of (D — 4) 
for the space-time dimension D tending to the physical value D = 4. As the 
result of a remarkable cancellation the third order pole in D — 4 disappears 
in the trajectory and the expansion starts with (D — 4)~ 2 . 
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1. INTRODUCTION 

Dynamics of semihard processes (i.e. processes with the c.m.s. energy yfs much 
larger than a typical virtuality Q) is one of the most interesting problems of modern 
high energy physics. There are a lot of papers devoted to the theoretical description 
of the processes (see, for example, Refs. [1,2] and references therein). Neverthe- 
less, the present situation here is not quite satisfactory, even in the framework of 
perturbative QCD which, presumably, can be applied to such processes |J. The 
BFKL equation || used for the calculation of the small x — — behaviour of cross- 
sections is responsible for the leading ln(l/x) terms only. To define a region of its 
applicability as well as to fix a scale of virtualities of the running coupling constant 
a s (Q 2 ) in the equation one has to know radiative corrections to the kernel of the 
equation. The crucial point in the programme of calculation of the next-to-leading 
corrections [|J is the gluon Reggeization in QCD. The Reggeization formed a basis 
of the original derivation of the BFKL equation M and was proved in the leading 
ln(l/x) approximation for elastic and inelastic processes 0. In this approximation 
the gluon trajectory 

j(t) = l + u(t) , 

u(t)=U W (t)+J 2 \t) + --- (1) 

can be written as an integral in the transverse momentum space: 
fll , . gH N r d D ~ 2 k ^ 9 

corresponding to a simple Feynman diagram of (D — 2)-dimensional field theory. 
Here and below N is the number of colours (N = 3 for QCD). The integral is 
divergent in the physical case D = 4. In fact one gets 

where 

D , g 2 N , 

e=y-2, f = - 6) . (4) 



This divergency is effectively cancelled in the kernel of the BFKL equation 

fc(q,Q ') = 2uj(t)S(q- q) +K, rea i(q, q ') , (5) 

9 ') = ^/(l + e) (g-V) 2 ' (6) 
which is not singular at D = 4 for integration with any smooth function of q ' . 

In the next-to-leading approximation Reggeization was checked and the cor- 
rection to the trajectory was presented in the form of two-loop integrals in the 
transverse momentum space [6,7]. However, it appears that we have to calculate 
the integrals because of increased infrared divergency of separate terms in the cor- 
rection to the kernel which diverge as (D — 4)~ 3 . So, contrary to what happened in 
the leading order, where the divergency cancellation was organized by a suitable re- 
arrangement of virtual and real contributions to the kernel without going to D ^ 4, 
in the correction one hardly can hope to do it. Therefore, we need to calculate the 
correction to the trajectory in an explicit form. In the present paper we solve this 
problem. 

2. REPRESENTATION OF THE QUARK AND GLUON 
CORRECTIONS 

The two-loop correction contains quark as well as gluon contributions. 

Let us first consider the quark contribution. Of course, here the case of light 
quarks is the most important one. Due to the cancellation of infrared terms in the 
equation masses of the light quarks should not affect the BFKL dynamics, so we 
can put them equal to zero. Therefore, according to Refs. [6,7] this contribution 
can be written as 

^ = I S^W jf dx{x{1 - X))W «* 2)e - 2( ^ )e ) • 

(7) 
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where n/ is the number of light flavours. Making use of the generalized Feynman 
parametrization 



r (i> 

at \i=l 



i=i 



n r i 



n 



\i=l 



*(i-£ 



.1 ; 



J][ / (iTj.v) 



a-i—l 



8=1 



E 

\i=l 



n 

ii=i 



(8) 



the integration in Eq. ([7p can be performed without any difficulty and leads to 



u>f\t)=g\q 2 ) 



2C 



4n f T 2 (2 + e) 



Ne T(4 + 2e) 
For the case e — > Eq. (§) gives 

^ 2) W = -^(g 2 ) 2e2n/1 



2 T 2 (l + e) 3 r(l-2e)r(l + e)r(l + 2e) 
7r(l + 2e) ~ 7r 2 (l-e) r(l + 3e) 



(9) 



3iV e 2 



(10) 



The gluon contribution is more complicated. It reads [0 



^a\t) = OfT {[^(1) + 2V(e) - ^(1 - e) - 2V(1 + 2e) 



1 



1 



1 + e 



(Jx - 2 JO - 2J 2 + 2 J 3 - J 2 



(11) 



(l + 2e) \e 2(3 + 2e) 
where ip(x) is the logarithmic derivative of the gamma function, ip (z) =d In T(z)/dz 
and the quantities Ii, J% are two-loop integrals in the transverse momentum space 
They can be written in the following form: 



9 2 N 



2(2n 



,£>-! 



2\6-D 



with 



and 



Ji 



ai = 1 



9 2 N ' 

D-l 



2(2< 



S D - 2 kid ( - D - 2 k2(q , 
Qi 2 (Qi - Q) 2 Q2 2 (Q2 -q) 2 

a2 = _ f^] 



rf( D - 2 ) gi rf( D - 2 )g 2 (g 2 ) 5 - D ; 
Qi 2 q2 2 (Q- Qi - Q2) 2 



(12) 



(13) 



(14) 
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with 

fe = 1 , „ = * = !»(£). (15) 

Let us note that because of the logarithmic factors (coming from the integration 
over longitudinal momentum components) the integrals Jj and Jj, with 2 > 2, do 
not correspond to usual Feynman diagrams of D — 2-dimensional field theory. 

The integral I\ is given by the square of the integral for the leading contribution 
to the trajectory (g) and reads 

Il = - 9 <(1^±A) 2 . (16) 

y V er ( 1 + 2e )/ 

The calculation of the integral J\ is not much more complicated than that of I\. 
It can be performed by the subsequent integration over q\ (using the Feynman 
parametrization to join the denominators q x and (q — qi — Q2) 2 ) and g*2- At the 
second step a non-integer power of (g*2 — <i) 2 appears in the denominator, so that 
we need to use the generalized parametrization (H). One easily obtains 

4 3 T(l - 2e)r 3 (l + e) 
J i-9 e 2 r 2(l- e )r(l + 3e) ' 1 ' 

The integrals J2 and J3 can be calculated in the same way with the help of the 
representation 

lna = — — a -1/ L = o . (18) 
av 

J 2 = J x [^{\ - e) + V>(2e) - V(l - 2e) - V(3e)] , (19) 
J 3 = Ji[^(l) + ^(e) - ^(1 - 2e) - ^(3e)] . (20) 

The calculation of the integral I2 is much more complicated and it appears that 

in the general case D 7^ 4 it can be given only in terms of infinite series. [] Instead 

1 We are thankful to A.V. Kotikov who paid our attention on papers || where I2 is expressed 
in terms of generalized hypergeometric series 



The result is 



in the interesting for physical applications case D — > 4, the integral is expressed 
in terms of the well-known values of the Riemann zeta function. Because of the 
importance of its contribution we present the calculation of the integral in the 
following section. The result is 



I - 9 

io — — 



-- + 2e^'(l) - 13eV'(l) 
e 



(21) 



Let us remind that 



V(1) = C(2) 



7T 



<(1) = -2C(3) 



-2.404 



(22) 



Expanding for the physical case e — > the integrals l\ and J1-J3, given respectively 
in Eqs. (0), flT7|), © and fl20|), and the coefficient of ii - 2Ji in Eq. (JJ) we find 

4^ 4 



Ii 



e* 2 
3^ 4 



J-2 



[1 - 2eV(l) - 2eV(l)] 



1 - 2eV(l) - 5eV (1) 
~ + e^(l) - f 6V'(1) 



J. 



s 4 r 2 



3 — 



-- + 4e^'(l) - 8eV'(l) 



(23) 

(24) 
(25) 
(26) 



and 



+ 2^(e) - V(l - e) - 2^(1 + 2e) + 



;i + 2e) Ve 2(3 + 2e)/ 

-i-t -(*«-§)«- G^f)'- 

Finally substituting the results fl2l|), (p3|) -( p7|) into Eq. (p|) , we obtain the gluon 
contribution to the two-loop correction to the gluon trajectory: 

1 



^ ) (t)=g\q 2 Y t - 9 



11 ( 1 1 / \ 67\ / ... , 22 . . . 404 . ., 
y +(2^(l)- y ) e+ (^(l)- y ^(l) + ^) e 2 



(2f 
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3. CALCULATION OF THE INTEGRAL J 2 

We introduce the Feynman parameters x± and x 2 to join the denominators 
depending on qi and q 2 correspondingly in Eq. (|T2"D, and the parameter z to join 
the results of these two parametrizations. Then we integrate over q 2 keeping the 
difference A = q\ — q 2 fixed and after applying the relation arrive at 

h = l o 1 dzJ o 1 dx 1 j\x 2 (--^j f(u,D)\ v=0 , (29) 

where 

r (5 - -) 

f(u,D) = -^(g 2 ) (6 - D+ ^(l-z) (30) 

7T 2 

r d D - 2 A 

X J (A 2 )^[z(A 2 + xi<f - 2xiA • q) + (1 - z)x 2 q 2 - (z(A - x x q) - (1 - z)x 2 q) 2 f~% 

The integration over A with the help of the generalized Feynman parametrization 
with parameter x and the subsequent integration over x by parts yield 

with 

Q(x) = zxi(l - xi) + (1 - z)x 2 (l - x 2 ) + z(l - z)(l - x)(x 1 - x 2 f . (32) 
It is convenient to present the derivative of f(u, D) with respect to v in the form 

-^f{v,D)\ v=Q = f x + f 2 + h , (33) 



where 



h = ~\ K*(l - z))f(0, D) = - F(6 2 D) { %^))1~d ln ^(l " z)) , (35) 
/3 = -F(6 - DM1 - z)f-l ^ dx HI - • (36) 



Let us denote the contributions of the functions fi to the integral I 2 in Eq. (|29| ) as 



Xj. It is easy to see that 



X 



1 — 




( 



id \ r 2 (i- g ) 

2de) t 



/ 1 = / 1 [V(l) + 2^(e)-V(l-e)-2^(2e)] . (37) 



The calculation of X 2 and X 3 is not so simple. We will perform it in the limit 
e — > 0. The idea underlying the calculation is to extract those regions of the 
integration range which lead to singular contributions for e — > and to put e = 
in the remaining region. Considering the integrand in Eq. ([12]) and our way of the 
Feynman parametrization it is easy to understand that the singular contributions 
can give only vicinities of points Xi = and Xi = 1, with % = 1,2. Taking into 
account the symmetry of the function Q(x) of Eq. ([32] ) under the simultaneous 
substitutions x% <-> 1 — X\ and x 2 1 — £2, we can restrict ourselves to the region 
Xi > 0, Xi + x 2 < 1 taking its contribution twice. It is convenient to introduce new 
variables defined through 



In these variables Q(x) becomes 

Q(x) = tR(x) , 

R(x) = zX{l - At) + (1 - z)(l - A)(l - (1 - X)t) + (1 - x)z(l - z){2\ - ift , (39) 



x\ = At , 



x 2 = (1 - A)t . 



(38) 



and we have 




A 



1 



(40) 




(41) 



where 




(42) 
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S (,,,)^r(2-2 e) /^.-7;^_|_. (43) 

The singular points are t — 0, A = and A = 1. Therefore, we divide the integration 
region in four pieces: 

a) 0<t<5 , < A < 1 ; b) 8 < t < 1 , < A < <5; 

c) 5 <t<l , < 1 - A < <5 ; d) 5 <t<l , 5 <\<l-5 . (44) 
Here we have introduced an intermediate parameter S assuming that 

e"« < 5 < e n (45) 

for any fixed power n. Due to the first of these inequalities one can put e = in the 
region d) and due to the second one we can neglect t, A and (1 — A) in comparison 
with 1 in the regions a), b) and c) correspondingly. Of course, 5 is cancelled in the 
whole expression for g(z) and g(x, z). 

In the region a) we have independently on x 

R = z\ + (1 — z)(l - A) (46) 

so that the integration in Eqs. ( [4*2]) and ( f4lf ) is quite simple. Performing a suitable 
expansion in e we find 



g^{z) = T{l-2e) ^ / " - " 2 



(^ 1 + (i-^ 1 ) + r^ ln ( r ^)] . ( 47 ) 

r(1 ~ 2e) ,5 2 V- 1 (z 2 - 1 + (1 - z) 2 *- 1 ) . (48) 



Notice that in g^ a \z) we need to keep terms which were neglected in g^ a \x, z). 
In the region b) R(x) takes the form 

R( x ) = A + r(s) , r(x) = (l-z)[l- t(l - z(l - x))} (49) 
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and considerably simplifies for x = 1. It is convenient to start with the integration 
over A without deriving with respect to x in g(x,z). After that we can put e = 
in the contributions containing 5 + r(x). The integration over t in the case of g(z) 
can be done without problems and gives 



=r(l-2e) 



} \V{At) 2e) l-z+6 { 5 2 



(50) 



As for g(x, z), using the approximate relation 



2e-i d x e e t 1 x e z{l + 2elni) d . , s^e-i (r-\\ 

1 dx~{r{ x )y- 2e = x{i- zy- 2e {i-tx) + i-z{i-x) Jv r ^ x >> ( ' 

the integration over t yields 

'1 — z\ 2t ~ x zx e 



g< b \x,z) = T(l-2e) 



+ 1(1 - zf^ In (I^) + j— ln(5 + s(l - z){\ - x)) 



(l-z(l-x)) 

1 (52) 

The calculation of g^ c \z) and g^ c \x,z) is immediate. In fact, because of the 
symmetry of the function R(x) in Eq. (^) under the replacement z ^ 1 — z 
and A <-» 1 — A, these two contributions can be obtained from Eqs. fl50D and (|52]) 
correspondingly by the substitution z <-> 1 — z. 

Finally we pass to the region d). After putting e = the integration over t 
becomes quite straightforward. The subsequent integration over A is very simple in 
the case of g(z) and gives 

* , n i £+*) + , n f Iz£±f) + 2 ,„ ( lr£) . (53) 

w 2 + 5 V 5 / 1-^ + 5 \ 8 J l-2z \ z J 
For the case of g(x, z) it is convenient to integrate without taking the derivative 
with respect to x. After integration over t one can represent the coefficient in the 
term with the logarithm as 

"^"-"'-^'-s Lu'-td-J (54) 



and integrating by parts over A arrives at 

Using the results given in Eqs. fl47D , (|50| ) and (|53|) we obtain with the required 
accuracy 

g(z) = g^( z )+g^(z)+g^(l-z)+g^(z) = T(l-2e)^± (z 2e ~ l + (1 - z) 2 *' 1 ) . 



Analogously, from Eqs. fl48|), (|52|) and ( |55"D for g(x,,z) we have 

g(x, z) = gW(x, z) + gO\x, z) + g®(x, 1 - z) + gW(x, z) 



(56) 



T(l - 2e) ((1 - ^) 2e -^ 2e + (1 - ^) 2e z 2e - 1 ) 

c e (l -x) 26 " 1 + - f31n(l-z) + -lnx) + 
xv 2 J 



(57) 



x V" 2 J 2xJ 
Substituting expressions (|56|) and (|5?|) into Eqs. (|40[) and fl4l|) respectively and 
performing the integration over x and £ we obtain 



J 2 = -2g 



r(l-2e) r 2 (2e)r(e)r(i- e ) 
p(i-c)r(4e) r(i) 

2^ 4 n 



(e)+ V(l-e) -2^(1)] 



- 2eip'{l) - 10eV(l) 



(5f 



and 



X 3 = -4^ 4 



r(l-2e)r(e)r(l + e) 
r 2 (l-e) r(l + 2e) 



g 4 n 



+ 2e^'(l) + 13eV(l) 



(59) 

we obtain for 



Using last two equations and Eq. (|37|) with given by Eq. ( 

J 2 = Xx + J 2 + X 3 the result (|T]). Q 

2 A.V. Kotikov informed us that he also obtained the same result from the generalized hyper- 
geometric series B 
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4. RENORMALIZATION AND DISCUSSION 

Up to now we worked with unrenormalized quantities. But renormalization is 
quite trivial. Since the trajectory itself must not be renormalized, we have only to 



use the renormalized coupling constant instead of the bare one g . In the MS 
scheme one has 



9 = 

where 



e 


1 + 


(t 


2n f \ 


9l' 








3 N J 


2e 



(60) 



2 ^AT(l-e) 
% = (47r)2+ e > ( 61 ) 

and for the gluon trajectory in the two-loop approximation we obtain 

2\ e of ^2 T /1 1 ^2 \ 



3 3 A^y \ 6 / 



^/_5 e+ /28_^ e2 y 



(62) 



3N \ 3 V 9 3 
The remarkable fact which occurred is the cancellation of the third order poles in 
e existing in separate contributions to the gluon correction ujf\t) ([□]). Possibly 
this fact indicates that the representation (|TT|) is not the best one. From the other 
hand a simpler representation for u^\t) means that there should exist a simple 
expression for the integral I2. In any case, the cancellation of the terms with e~ 3 
is very important for the absence of infrared divergences in the corrections to the 
BFKL equation. As the result of this cancellation the gluon and quark contributions 
to uj^{t) have similar infrared behaviour. Moreover, the coefficient of the leading 
singularity in e is proportional to the coefficient of the one-loop (3 function. This 
means that infrared divergences are strongly correlated with ultraviolet ones. The 
correlation is unique in the sense that it provides the independence of singular 
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contributions to u(t) on (f. Indeed, expanding Eq. 



u(t) = 
'67 



+ 



9 



-si 



7T 

y 
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y 



2tv 
3 N 



we have 
/ 1 



-In* 



V' 



ion/ 

9 N 



2 In 



404 > , , 56 n f 
17 + 2C(3)+ 27i 



(63) 



Eq. ([63]) exhibits explicitly all singularities of the trajectory in the two-loop ap- 
proximation and gives its finite part in the limit e — > 0. Of course, the singularities 
should cancel putting all the corrections into the BFKL equation. We hope to show 
this in a subsequent paper. 
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